Deceleration mechanism of Fermi acceleration in a time-dependent stadium billiard 
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The dynamics of a time-dependent stadium-like billiard are studied by a four dimensional nonlin- 
ear mapping. We have shown that even without any dissipation, the particle experiences a decrease 
on its velocity Such condition is related with a critical resonance velocity where if the initial ve- 
locity has a higher value than the resonant one, we can observe Fermi acceleration, however, if the 
initial velocity has a initial value smaller than the critical one, the particle is temporarily trapped 
surrounding the stability islands, in a stickiness regime. We believe that this sticky orbits can act 
as deceleration mechanism for Fermi Acceleration. 
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Problems involving particles colliding (with or without 
interaction among themselves) into a closed domain, are 
well known in the literature as billiard problems. The 
laws of reflection of the point particle (billiard ball) are 
the same for the light's specular reflection, i. e. while the 
normal component changes its sign on the collision point, 
the tangent one is preserved and the particle is reflected 
elastically [l[ . The study of billiard problems had begun 
with Birkhoff [2], but it was Sinai [3| that gave us the 
necessary mathematical formulation for a new rigorous 
analysis. Nowadays, billiard problems can be found in 
several fields of physics and other sciences, like optics 
0], microwaves [Bj, quantum dots ultra-cold atoms 
trapped in a laser potential 0, among others. 

If the billiard boundaries are time-dependent, we may 
found an interesting phenomenon called Fermi Accelera- 
tion (FA). Introduced for the first time in 1949 by Italian 
physicist Enrico Fermi [9[, FA is basically the unlimited 
energy growth for a classical particle suffering elastic col- 
lisions with a time-dependent boundary. For one dimen- 
sional billiards we have some examples of this unlimited 
growth of energy [Iol4l4j. In two dimensional billiards, 
in which the dynamics are more complicated, this phe- 
nomenon has already been studied in different models 
[l5l-[2Cij. One of the main questions about FA in per- 
turbed systems is whether the energy of the particle can 
grow to infinity. The answer is far away from trivial, and 
it depends on the kind of perturbation on the boundary 
as well the boundaries geometry. The Loskutov- Ryabov- 
Akinshin (LRA) conjecture [21[ tell us that if a billiard 
with fixed boundary presents chaotic behavior, Fermi ac- 
celeration is exhibited, when perturbation on the bound- 
ary is introduced. However, F.Lenz et. al [22j, shown 
that even for an elliptical billiard, which has regular dy- 
namics in the static boundary case, exhibits FA when 
time-dependence is introduced. In this particular case, 
FA production mechanism is due orbits that crosses the 
stochastic layer, changing their dynamics from librator 
to rotator, or vice- versa. And recently, Leonel and Buni- 
movich [23|, extended the LRA conjecture, by showing 
that only the presence of heteroclinc orbits in the static 
phase space, is sufficient condition for the system present 



FA when time-dependence is introduced. 

There are some methods to suppress FA, like intro- 
ducing dissipation via inelastic collisions [3, EH, , or 
through a drag force [20]. But, in this letter, we found 
that orbits in a stickiness regime [24], characterized 
by a resonance, acts as a deceleration mechanism of the 
particle average velocity. This stickiness phenomenon has 
already been studied in some billiards 113, HI E3 with 
different approaches, but in this letter we intent to show 
how this stickiness, can be used as a deceleration mech- 
anism for FA, even with no dissipation introduced at all. 

In this letter, a stadium-like billiard is studied con- 
sidering the dynamics proposed by Loskutov et. al. [2~0 |, 
which describes the billiard dynamics by a nonlinear map 
and considered the fixed boundary approximation, i. e., 
the boundary is considered as fixed, but when the parti- 
cle collides with it, they exchange momentum, as if the 
boundary were moving through a periodic time pertur- 
bation. A resonance critical value for the particle's veloc- 
ity, depending on the geometrical control parameters of 
the billiard is found by linearizing the unperturbed map. 
Our results show that depending on the initial velocity of 
the particle we can found FA for initial velocities above 
this critical value. For values of initial velocities lower 
than the critical one, we observe a decreasing on velocity 
even with no dissipation. This decreasing is related by 
some orbits that rotates around the fixed points due a 
resonance on the boundaries perturbation, in a regime 
of sticky orbits. As a consequence, the particle average 
velocity experiences a deceleration. 

The model under study consists of a classical particle 
suffering elastic collisions inside a stadium-like billiard 
with focusing components that are periodically time- 
dependent according to B(t) = Bo cos(u>£), where w is 
the frequency of oscillation and Bo is the amplitude of 
oscillation. The geometrical control parameters are de- 
scribed in FigHJ Suppose that focusing components, 
which are symmetric about the vertical billiard axis, had 
a radius R with the angle measure 2^>. Looking at the 
FigJTJ one can geometrically obtain the billiard parame- 
ters as R = (a 2 + 4b 2 ) /8b and $ = arcsin(a/2i?). 

The initial angular conditions are (a, <p). Just for no- 
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tation, all the variables with a star index, are measured 
just before the collision point. We assume that V n is the 
particle velocity and t n is the time of the n th collision, 
and of course, at initial time to = 0, the particle belongs 
to the focusing component and the velocity vector di- 
rects towards to the billiard table. In order to describes 
the dynamics of this billiard, we should consider two dif- 
ferent cases. In the first one, the particle collides with 
a boundary component and the subsequent collision is 
with the same focusing component. The second case, 
the particle collides with the focusing component, and 
in the next collision, the particle hits opposite boundary 
component, so in this case we use the unfolding method 
[H, EH to describes its dynamics. For both cases, the re- 
currence relation for the velocity and the angle a n are 
the same. Making a vectorial arrangement for V n +i and 
V n with the angles a n and a* and applying the cosine 
law for the velocity vectors V n , V n +\ and B ni according 
to FigJT]we found: 

a n = arcsin \-^— sin«) J . (1) 
V Vn+l J 



V n+1 = Vn + 4B n 2 + W n B n cos«) , (2) 

where B n is the boundary velocity vector in the time t n . 

Let us now introduce the dynamical equations for the 
first case, i.e., the successive collisions with the same fo- 
cusing component. This kind of collisions will only hap- 
pen when <l> >| (p n +i |. According to the figure [2] and the 
specular reflection law, we can obtain these relations: 

a n+l — a n 

^n+i = ^ n + 7T - 2a n (mod 2tt) . 




(3) 



FIG. 1: (color online). The unfolding method, the geometrical 
control parameters and a trajectory with indirect collisions. 

Now, we should consider the case where the particle 
has its next collision with the opposite focusing compo- 
nent. For this case, it is necessary that $ <| <^ n +i |, and, 
we will introduce some more dynamical variables, as i/j, 



that is the angle between the particle trajectory and the 
vertical line at the collision point, and x n , that is the 
projection of the particle position under the horizontal 
axis. 

Looking at the figure [2] we can geometrically obtain the 
value of the angle ip n = a n — (p n . We also can see that x n 
is the sum of the line segments AB + BC + CD. Taking 
into account the value of ip n , and, after some geometrical 
algebra, we obtain x n = CQ8 f^ n) [sin(a n ) + sin(^> - ^ n )). 
The recurrence relation between x n and x n+ i is given 
by the unfolding method, described in FigJTJ as Xrr,A_~\ — 
x n + /tan(VVi)- 

Now let us find the equations for the angular dynami- 
cal variables and time. If we invert the particle motion, 
i.e. consider the reverse direction of the billiard parti- 
cle, then the expression that furnishes us the value of x n 
is also inverted and the angle a n: became a*. Resolv- 
ing it with respect to a*, taking into account that this 
angle is changed in the opposite direction than a n , and 
the angle cp n should have the reversed sign, we now can 
get the value of the incident angle a* , that will become 
Q^+i wnen we re-inverted the particle motion. The val- 
ues of ifn+i and the time t n+ i can be obtained by easy 
geometrical considerations on FigJTJ Thus, we obtain the 
mapping the case where collisions with the opposite fo- 
cusing component happen. 

= arcsin[sin(^ n + $) - ^n+i cos(^ n )/^] • 

R[cos((p n ) + cos(y? n+ i) - 2 cos(^)] + / 

— + , . v 

KcOS(^ n ) 

(4) 




FIG. 2: (color online). Some dynamical variables and a tra- 
jectory with successive collisions. 

When we linearize the unperturbed map [2l|, [28| , we 
may obtain, according the action-angle variables, a rota- 
tion number as a = arccos(l — 8bl /a 2 cos 2 (^)), where 
ip n = arctan(ma//) is the fixed point, where m > 1, is the 
number of mirrored stadiums in the unfolding method. 
Considering a trajectory where the particle moves around 
some stable fixed point, the time between two sequential 
collisions is given by r ~ 1/ cos(^*)V^. Thus, the ro- 
tation period is T rot = 27rr/<j. If the rotation period is 
equal to the period of external perturbation T ext = 2tt/u, 
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we observe a resonance between rotation and boundary 
oscillations. This resonant velocity is given by: 



r cos(^*) arccos(l — Sbl/(a cos(^*))) 

It is always good to remember that, this resonant phe- 
nomenon only occurs when the defocusing mechanism 
[29| in this billiard is not working, according to the ex- 
pression 1/2R w 4bl/a 2 > 1. 

If the initial velocity of the particle is higher than the 
resonant one, we observe typical FA behavior when we 
study the evolution of V as function of the number of 
collisions. Figure [3j shows some evolutions of V for 5000 
different initial conditions (a, (/?), evaluated up to 10 9 col- 
lisions, for Vq > V r , where j3 ~ 0.5 is the growth expo- 
nent. However, not all initial conditions are accelerated. 
In the dark part of Figj3]^b) are the initial conditions who 
suffered the effects of FA. On the other hand, the initial 
conditions that are inside the islands, do not experience 
this unlimited energy growth. We evaluated each pair 
(a,(p) up to 10 8 . If after this time, the particle has its 
velocity increased at least one order of magnitude, we 
considered that the particle suffered FA. As far we can 
study, the initial conditions who not suffered FA, are in 
the regular region of the phase space shown in FigH^a). 
It is important to clarify, that the control parameters 
used in Figj3j will be the same used in all figures of this 
letter. 




(a) n (b) 



FIG. 3: (color online), (a) Evolution of some curves of V 
as function of the number of collisions for Vb > V r . The 
control parameters used were a — 0.5, b — 0.01, / — 1, Bo — 
0.01, which give us V r = 1.2. (b) Grade of 500 x 500 initial 
conditions (a, ip G [0, tt/2]) for Vb = 5. 

However, if the initial velocity has a smaller value 
than the resonant critical one given by Eq.(j5j), the par- 
ticle experiences a resonance, and stays for limited time, 
surrounding the stability regions. In this time interval, 
the particle can penetrate into the neighborhood of fixed 
points, as a result, the whole region of the phase space 
became accessible [28], as shown FigHJwhere some phase 
spaces are shown for different values of Vb- It is clearly 
to see, that if Vb > V r , the particle does not experience 
the resonance and the phase space does not become all 
accessible, differently in what happen if Vb < V r . 

This resonance can appear in systems with mixed 
phase space [29[ and it causes a stickiness phenomenon 




FIG. 4: (color online) Phase Space for 25 different initial 
conditions iterated up to 10 6 collisions. In (a) Vb = 5, (b) 
V = 1.3, (c) V = 1.2 and (d) V = 0.5. 

[zil [25j . where some orbits can be trapped, surround- 
ing stability islands for a limited time interval. We do 
believe, that this orbits in stickiness regime, are the re- 
sponsible for the deceleration mechanism of the velocity. 
Looking at Figj5^a,b) where a single initial condition is 
iterated up to 10 5 collisions, one can see, that after the 
orbit experiences a initial stickiness regime (colored parts 
in both Figs|5fa,b)), the velocity has a decrease. 




FIG. 5: (color online) In (a) we show an orbit with an initial 
stickiness and trajectories penetrating through the stability 
islands for Vb = 0.5; and in (b) Velocity evolution as function 
of n, for the same initial condition of (a), where the colored 
parts, corresponds to the same stickiness behavior shown in 
(a). 

This kind of stickiness influence is interesting and need 
more attention, since this kind of dynamics is considered 
an open problem [30]. We can see better this effect in 
Figj6]^a,b,c) where we show the dynamics evolution in a 
different phase space where the coordinates are velocity 
and the time taken mod(27r) for 50 initial conditions iter- 
ated up to 10 5 . We can see some darker lines in the region 
of the initial velocity, indicating the orbits in stickiness 
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FIG. 6: (color online) Phase Space in the coordinates V and 
t mod{2n) for 25 different initial conditions iterated up to 
10 6 collisions. In (a) V = 0.6, (b) V = 0.5, (c) V = 0.4 and 
(d) a resonance zoom-in with Vb = 0.9. 



regime in the range of the initial velocity. A zoom-in 
over these dark curves is shown in FigE^d), where we 
are able to see the rotations, during this initial stickiness 
behavior. 

When we consider the average velocity, we will found 
similar behavior as shown FigEl The V curves experi- 



ences a constant velocity initial dynamics, where the or- 
bits are experiencing this stickiness behavior confirmed 
by FigsE] and [6j and then they start decreasing and 
reaches a constant plateau, remaining in this dynam- 
ics for long periods (in our numerical simulations up to 
3 x 10 6 ). 




FIG. 7: (color online) Behavior of the average velocity over an 
ensemble of 5000 different initial conditions, for some different 
values of initial velocity. The velocity decreasing is evident. 

To conclude, we have presented that orbits in stickiness 
regime can act as a deceleration mechanism for FA in a 
time-dependent stadium billiard, even with no dissipation 
at all. At principle, we expect that this mechanism can 
be extended to other time-dependent billiards. 

ALPL and AL acknowledges FAPESP for financial 
support and E.D.L. acknowledges CNPq, FAPESP, 
CAPES and FUNDUNESP. 



[1] N. Chernov and R. Markarian, Chaotic Billiards. (Amer- 
ican Mathematical Society, Vol. 127, 2006). 

[2] G. D. Birkhoff, Dynamical Systems. (Amer. Math. Soc. 
Colloquium Publ. 9. Providence: American Mathemati- 
cal Society. 1927). 

[3] Y. G. Sinai, Russ. Math. Surveys 25, 137 (1970). 

[4] E. D. Leonel, Phys. Rev. Lett. 98, 114102 (2007). 

[5] J. Stein and H. J. Stokmann, Phys. Rev. Lett. 68, 2867 
(1992). 

[6] C. M. Marcus et a/., Phys. Rev. Lett. 69, 506 (1992). 
[7] N. Friedman et a/., Phys. Rev. Lett. 86, 1518 (2001). 
[8] M. F. Andersen et a/., Phys. Rev. Lett. 97, 104102 
(2006). 

[9] E. Fermi. Phys. Rev. 75, 1169 (1949). 
[10] A. J. Lichtenberg, M. A. Lieberman, Regular and Chaotic 
Dynamics, (Appl. Math. Sci. 38), Springer - Verlag, New 
York 1992. 

[11] E. D. Leonel, J. K. Leal da Silva and P. V. E. McClintock, 

Phys. Rev. Lett. 93, 014101 (2004) 
[12] E. D. Leonel and P. V. E. McClintock, J. Phys. A 38, 

823 (2005). 

[13] E. D. Leonel and A. L. P. Livorati, Physica. A 387, 1155 
(2008). 

[14] A. L. P. Livorati, D. G. Ladeira and E. D. Leonel, Phys. 

Rev. E 78,056205 (2008). 
[15] D. F. M. Oliveira and E. D. Leonel, Physica A 389, 1009 



(2010). 

[16] E. D. Leonel, D. F. M. Oliveira and A. Loskutov, Chaos 
19, 033142 (2009). 

and Numer Simulat. 15, 1092 (2010). 
[17] F. Lenz et a/., New J. Phys, 11, 083035 (2009). 
[18] J. A. de Oliveira, R. A. Bizao and E. D. Leonel, Phys. 

Rev. E, 81, (2010). 
[19] A. K. Karlis et al, Phys. Rev. Lett. 97, 194102 (2007). 
[20] E. D. Leonel and L. Bunimovich, Phys. Rev. E, 81, 

(2010). 

[21] A. Loskutov, A. B. Ryabov and L. G. Akinshin, J. Phys. 

A, 33, 7973 (2000). 
[22] F. Lenz, F. K. Diakonos and P. Schemelcher, Phys. Rev. 

Lett. 100, 014103 (2008). 
[23] E. D. Leonel and L. Bunimovich, Phys. Rev. Lett. 104, 

224101 (2010). 

[24] G. M. Zaslasvsky, Hamiltonian Chaos and Fractional Dy- 
namics, Oxford University Press, New York (2008). 

[25] G. M. Zaslasvsky, Phys Rep., 371, 461 (2002). 

[26] E. G. Altmann, A. E. Motter and H. Krantz, Phys. Rev. 
E, 73, 026207 (2006). 

[27] E. G. Altmann, A. E. Motter and H. Krantz, Chaos, 15, 
033105 (2005). 

[28] A. B. Ryabov and A. Loskutov, J. Phys. A, 43, (2010). 
[29] L. Bunimovich, Commun. Math. Phys. 65, 295 (1979). 
[30] L. A. Bunimovich, Nonlinearity, 21, 13 (2008). 



